Perturbations of Schwarzschild Black Holes in Chern-Simons Modified Gravity 
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We study perturbations of a Schwarzschild black hole in Chern-Simons modified gravity. We begin 
by showing that Birkhoff 's theorem holds for a wide family of Chern-Simons coupling functions, a 
scalar field present in the theory that controls the strength of the Chern-Simons correction to the 
Einstein-Hilbert action. After decomposing the perturbations in spherical harmonics, we study the 
linearized modified field equations and find that axial and polar modes are coupled, in contrast 
to general relativity. The divergence of the modified equations leads to the Pontryagin constraint, 
which forces the vanishing of the Cunningham-Price- Moncrief master function associated with axial 
modes. We analyze the structure of these equations and find that the appearance of the Pontryagin 
constraint yields an overconstrained system that does not allow for generic black hole oscillations. 
We illustrate this situation by studying the case characterized by a canonical choice of the coupling 
function and pure-parity perturbative modes. We end with a discussion of how to extend Chern- 
Simons modified gravity to bypass the Pontryagin constraint and the suppression of perturbations. 

PACS numbers: 04.50.Kd,04.25.-g,04.25.Nk,04.30.-w 



I. INTRODUCTION 

Extensions of general relativity (GR) are inherently 
interesting because they hold the promise to address un- 
resolved problems in cosmology and astrophysics. One 
such extension is Chern-Simons (CS) modified gravity 
which has recently been used to propose an explanation 
to the cosmic baryon asymmetry Q and polarization in 
the cosmic microwave background (CMB) Q. This ex- 
tension has also attracted some recent interest because it 
might allow for a test of a wide class of string theories [H • 

CS modified gravity introduces a well-motivated cor- 
rection to the Einstein-Hilbert action: the product of a 
parity-violating, Pontryagin term (59j and a scalar field 
or CS coupling function that controls the strength of 
the correction. Such a modification has deep roots in 
the standard model, since chiral, gauge and gravitational 
anomalies possess Pontryagin-like structures. CS modi- 
fied gravity is also motivated by string theory and loop 
quantum gravity [f|. In the former, it arises through the 
Green-Schwarz mechanism [6|, as an anomaly-canceling 
term. In fact, the CS term is a requirement of all 
model-independent extension of 4-dimensional compact- 
ifications of string theory Q • 

Some theoretical aspects of CS modified gravity have 
recently been investigated. Cosmological studies have 
been carried out in [3|, |8|, |9| in connection with the Cos- 
mic Microwave Background radiation and in 0, fTol ] in 
the context of leptogenesis. Gravitational wave solutions 
have also been studied in 0,13,110, where they were 
found to possess amplitude birefringence, possibly lead- 
ing to a test of the theory Q. Weak- field solutions for 



spinning objects in CS modified gravity have been stud- 
ied in [12 . Il3| ], leading to a prediction of gyromagnetic 
precession that differs from GR, and later to a constraint 
on the magnitude of the CS correction [l^|- Further dis- 
cussion on these issues and others related can be found 

in iMP M EE M E3, M M IE HE IE IE Hi IE 

12a . [271 l28j and references therein. 

Solutions to the CS modified field equations that rep- 
resent interesting physical configurations have also at- 
tracted attention. Apart from first formalizing the the- 
ory, Jackiw and Pi [l[ also showed that the Schwarzschild 
and gravitational plane-wave spacetimes remain solu- 
tions in CS modified gravity. The Reissner-Nordstrom 
and Friedman- Robertson- Walker metrics were also found 
to persist in the theory (24j. On the perturbative 
front, H HH found a weak-field solution for orbiting, 
spinning objects, while |25[ showed that the weak-field 
limit of the Kerr metric remained a CS solution to first 
order in the metric perturbation for a specific choice non- 
canonical of the CS coupling. Exact solutions that rep- 
resent spinning objects in CS modified gravity have so 
far only been studied in _29] , where it was suggested that 
these exist provided the gravitomagnetic sector of the 
metric is non-vanishing or that stationarity is broken. 

In this paper, we study perturbations of a 
Schwarzschild black hole (BH) in CS modified grav- 
ity We begin by showing that Birkhoff 's theorem - 
the statement that the Schwarzschild solution is the only 
vacuum, spherically-symmetric solution of the theory- 
persists in the modified theory for a wide class of cou- 
pling functions. After decomposing the metric pertur- 
bations in spherical harmonics, we find equations that 
govern the behaviour of each harmonic. In GR, the di- 
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vergence of the field equations (both at the perturbative 
and non-perturbative levels) implies energy-momentum 
conservation and, in the absence of matter, it becomes 
an identity that ensures diffcomorphism invariance. In 
CS gravity, however, the divergence of the field equations 
(again at the perturbative and non-perturbative levels) 
leads to a constraint, the so-called Pontryagin constraint, 
which ensures that the theory remains diffeomorphic in- 
variant and that the Strong Equivalence Principle holds. 
This constraint imposes a restriction on the class of met- 
rics that can be solutions of the theory, specially in vac- 
uum [2^] . At the perturbative level, the Pontryagin con- 
straint is automatically satisfied for a flat background, 
but we show that for a BH background it forces the 
Cunningham-Price-Moncrief (CPM) master function, as- 
sociated with axial perturbations (perturbations of odd 
parity), to vanish. Such a strong restriction on the pos- 
sible perturbative modes a BH is allowed to possess is a 
distinctive feature of CS gravity. 

After investigating these preliminary issues, we con- 
centrate on the study of general perturbations of a 
Schwarzschild spacetime [6l|. We find that the CS mod- 
ification introduces new terms into the field equations 
that results in a mixing of polar and axial parity modes. 
Therefore, in CS modified gravity modes with different 
parity do not decouple, as is the case in GR. Moreover, 
the new terms contain third-order derivatives that change 
the basic structure of the partial differential equations 
that describe the behavior of the metric perturbations. 
In particular, the Pontryagin constraint constitutes a new 
equation for the metric perturbations, which in general 
leads to an overdetermined system. A priori, it remains 
unclear whether this system has non-trivial solutions, 
i.e. whether the entire set of field equations is compati- 
ble. This paper shows that at least a wide class of BH 
oscillation modes are forbidden in the modified theory 
for a large family of CS coupling functions. More specif- 
ically, we show that it is not possible to have either pure 
axial or pure polar oscillations of a BH in CS gravity for 
a wide class of coupling functions. 

We end with a discussion of possible routes to extend 
CS modified gravity such that generic BH oscillations are 
allowed. A promising route is to allow the CS coupling 
function to be a dynamical quantity. In this case, the 
Pontryagin constraint is balanced by the equation of mo- 
tion of the scalar field, thus preventing the CPM function 
to identically vanish. Moreover, the field equations are 
modified by the introduction of a stress-energy tensor for 
the CS scalar field and by terms describing perturbations 
of the scalar field. The extended set of field equations, 
which has a certain ambiguity encoded in the potential 
of the CS scalar field, is no longer overdetermined, and 
hence, it provides a suitable framework to study the mod- 
ified dynamics of BH oscillations. 

This paper is divided as follows: Sec. |TT] describes 
the basics of CS modified gravity; Sec. |HT] establishes 
Birkhoff 's theorem for a family of CS coupling functions; 
Sec. IIVI begins by describing the basics of Schwarzschild 



BH perturbation theory, and then discusses the conse- 
quences of the linearized Pontryagin constraint and the 
structure of the modified field equations; Sec. |V] inves- 
tigates perturbations of canonical CS gravity, including 
the impossibility of having either purely polar or purely 
axial perturbative modes; Sec. I VII points to directions in 
which the theory could be extended to bypass the Pon- 
tryagin constraint and allow for generic oscillations of a 
Schwarzschild BH; Sec. IVIII concludes and points to fu- 
ture work. 

We use the following conventions throughout this work. 
Greek letters and a semicolon are used to denote in- 
dices and covariant differentiation respectively on the 4- 
dimensional spacetime. In some cases, we denote covari- 
ant differentiation with respect to a Schwarzschild back- 
ground metric by V„. Partial differentiation of a quan- 
tity Q with respect to the coordinate x M is denoted as 
d xtt Q or Q Symmetrization and antisymmetrization 
is denoted with parenthesis and square brackets around 
the indices respectively, such as A, b -> := [A ab + A ba ]/2 
and Aj ab j := [A ab — A ba ]/2. We use the metric signature 
(— , +, +, +) and geometric units in which G = c = 1. 



II. INTRODUCTION TO CHERN-SIMONS 
MODIFIED GRAVITY 

In this section we describe the basics of CS modified 
gravity as it was introduced by Jackiw and Pi [l| . (more 
details can be found in P, [H, and references therein). 
In this paper, we shall be mainly concerned with this for- 
mulation [lj, but we will also discuss possible extensions 
later on. 

The CS extension of the GR action is given by 



&cs '■— 



d A x, 



R + Anat - -9*RR 



(1) 



where K = l/(167r), g is the determinant of the space- 
time metric g„„, R stands for the Ricci scalar and £ mat 
is the Lagrangian density of the matter fields. We rec- 
ognize the first two terms as the Einstein-Hilbert action 
in the presence of matter, while the last term is the CS 
modification, defined via 

*R R := R a0jS *R a ^ s = -R a ^ lS r] al3 ^ v R lS ^ v , (2) 

where the asterisk denotes the dual tensor, constructed 
using the 4-dimensional Levi-Civita tensor r} a t 3 i lv [62| . 
The strength of the CS correction is controlled by the 
scalar field 6, which we shall refer to as the CS cou- 
pling function or CS scalar field. The so-called canonical 
choice of 9 corresponds to |l[ 



0can ~ [t/fi, 0, 0, 0] , 

where \x is a dimensional parameter. 



(3) 
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As usual, we obtain the (modified) field equations by 
varying the action with respect to the metric, yielding 

G p.v + c p.v = 87rT^ at , (4) 

where T™ at is the matter stress-energy tensor. We shall 
refer to the quantity C ^ as the C-tensor, defined as 

c^^cjn + cjn, (5) 

where 

c# O-^^Ru)^, (6) 

The spatial sector of the C-tensor reduces to the 3- 
dimensional Cotton tensor in some symmetric cases [63[ . 
For this reason, the quantity y := 6 is sometimes re- 
ferred to as the embedding coordinate, since it embeds 
the 3-dimensional CS theory into a 4-dimensional space- 
time. 

Diffeomorphism invariance is preserved provided an 
additional equation is satisfied. This equation can be 
obtained by computing the covariant divergence of the 
modified field equations 

C^ = ^v„*RR = &nT™ t >» (8) 

In vacuum, T'™ at = and thus we are left with the so- 
called Pontryagin constraint 

*RR = 0, (9) 

which is an additional equation that the metric tensor 
has to satisfy. 

III. BIRKHOFF'S THEOREM IN 
CHERN-SIMONS MODIFIED GRAVITY 

Birkhoff's theorem states that the most general, spher- 
ically symmetric solution to the vacuum Einstein equa- 
tions is the Schwarzschild metric. In CS modified grav- 
ity, this is not necessarily the case because the C-tensor 
modifies the field equations. In [l[ it was shown that the 
Schwarzschild solution persists in CS modified gravity, 
but this does not necessarily imply that Birkhoff's theo- 
rem holds. Let us then study Birkhoff's theorem in CS 
modified gravity in more detail. 

The line element of a general, spherically-symmetric 
spacetime can be written as a warped product of two 
2-dimensional metrics [30l l3lj : a Lorentzian one, g AB 
(A, B, . . . , H = t,r), and the unit two-sphere metric, £l ab 
(a,b, . . . , h — 9, (/>). This line element takes the following 
2 + 2 form: 

g dx»dx v = g AB (x c )dx A dx B + r 2 {x A )VL ab {x c )dx a dx h , 

(10) 



where the warped factor is the square of the area radial 
coordinate r. Covariant differentiation on the Lorentzian 
manifold is denoted by a bar while on the 2-sphere is 
denoted by a colon. The case of the Schwarzschild metric 
in Schwarzschild coordinates in given by: 

g AB dx A dx B = -fdt 2 + r 1 dr 2 , f = l-™, (11) 

where M is the BH mass. 

The Cotton tensor associated with this metric vanishes 
identically if the scalar field has the following form [&4l ]: 

9 = 9{x A )+r{x A )<d{x a ). (12) 

This functional form is invariant under coordinate 
changes that leave the 2 + 2 structure of the metric 
[Eq. (jlOp ] invariant, i.e. coordinate transformations: 
x A — ► x A — x A (x B ), x a — ► x a = x a (x b ), and 
transformations on the unit two-sphere. Moreover, one 
can show that the two pieces of the C-tensor shown in 
Eqs. © and ([7]) vanish independently for the family of 
scalar fields given in Eq. (TT2")) . For most of this paper, we 
will consider the particular case in which the scalar field 
also respects the spherical symmetry of the background: 

9 = 9{x A ). (13) 

Let us comment some more on this result and place it 
in context. The CS scalar field 9 is usually assumed to 
depend only on cosmic time, = 9(t). This assumption 
has its roots in the work of Jackiw and Pi [l[ , who fur- 
ther imposed that 9 = d9 j dt be constant, such that time- 
translation symmetry and space-time reparameterization 
of the spatial variables be preserved. With these assump- 
tions, they showed that CS gravity can be interpreted 
as a 4-dimensional generalization of 3-dimensional Cot- 
ton gravity, and that the Schwarzschild and Friedman- 
Robertson- Walker solutions persist in the modified the- 
ory. Later on, Smith, et. al. [3] argued that the CS 
scalar field could represent some quintessence field that 
enforces the arrow of time associated with cosmic expan- 
sion. In this paper we have considered the most general 
spherically-symmetric spacetime and we have written its 
line element in the 2 + 2 form that follows from its warped 
geometric structure. This form of the metric is invariant 
under general coordinate transformations in the two 2- 
dimensional manifolds, associated with this 2 + 2 struc- 
ture. Therefore, it is not too surprising that if Birkhoff's 
theorem is satisfied for the canonical choice of CS scalar 
field, it is also satisfied for fields that depend arbitrarily 
on t and r. 



IV. PERTURBATIONS OF A 
SCHWARZSCHILD BLACK HOLE 

A. Basics 

We begin with a short summary of the basics of metric 
perturbations about a Schwarzschild background. The 
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spacetime metric of a perturbed BH can be written in 
the form 



"flU 1 



(14) 



where g is the background Schwarzschild metric (ge- 
ometric objects associated with it will have an overbar) 
and h is a generic metric perturbation. Thanks to the 
spherical symmetry of the background, we can expand 
the metric perturbations in (tensor) spherical harmon- 
ics. In this way, we can separate the angular dependence 
in the perturbative equations, which yields a much more 
simple system of equations: a system of 1+1 partial dif- 
ferential equations (PDEs) in time and in the radial area 
coordinate r. In addition, we can distinguish between 
harmonics with polar and axial parity [651 ] . also called 
even and odd parity modes respectively. In GR, the per- 
turbative field equations decouple for modes of different 
parity, but this may not be the case for alternative the- 
ories. 

Let us then split the metric perturbation /i into polar 

and axial perturbations, h^ v = h^ v + /iJL, and each of 
these into (tensor) spherical harmonics via 



(15) 



where 
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• a,tm 



ulm -yl 
n AB 1 



t.m 



h 1 ^ s e a m 



. Trim cim 



Pa Y a 



(16) 



,(17) 



and where asterisks denote components given by sym- 
metry. The quantity Y im refers to the standard scalar 
spherical harmonics [sec 32] for conventions] , while Y^ m 
and S^ m are polar and axial, vector spherical harmonics, 
defined only for t > 1 via 



Y 



Im 



Y 

1 -n 



S. 



Va 1 b 



(18) 



Similarly, Y£™ and Z^™ are polar, and S^ 1 axial, sym- 
metric tensor spherical harmonics, defined only for £ > 2 
via 



\^£m \rtmc\ 

1 ab = 1 "ab > 



J ab 



Y. 



ab 



s, 



ab 



qIttl 



fclly^Q ob ,(i9) 



(20) 



The sign convention for the Levi-Civita tensor of the 2- 
sphere is: r\ s — sin 6. All metric perturbations, scalar 
{^ab)i vectorial {p ™ and q^ 1 ), and tensorial (K 
G lm , and q^™ 1 ), are functions of t and r only. 



In GR, the Einstein equations can be decoupled in 
terms of complex master functions that obey wave- 
like master equations. Once the master functions are 
constructed we can recover all remaining metric per- 
turbations from them. For axial modes, we can use 
the Cunningham-Price-Moncrief (CPM) master func- 
tion [33[, defined by 



A, 



htt - hi™ + -h 



' r,l 



t,r 



(21) 



whereas for polar modes we can use the Zerilli-Moncrief 
(ZM) master function 0, ^ 



A/™ — 

1 + A, 



{K em + (1 + \ t )G 



I'm 



J_ 

A, 



ffC - rK* 



-(l + \/)Pr 



t in 



,(22) 



where X g = (£+ 2)(£ - l)/2 and A e = X e + 3M/r. These 
two complex master functions are gauge invariant. In or- 
der to simplify our analysis we shall here fix the gauge by 
setting the following metric perturbations to zero (Regge- 
Wheeler gauge): 



jjtm = o j 

Glm ~&m 
— Pt 



= 0. 



(23) 
(24) 



The master equations for the master functions, in the 
case of perturbations without matter sources, have the 
following wave-like structure: 



v 



■/Ax 



ivn 

Cpm/zm 



= o, 



(25) 



coordinate r+ 

, r /Ax 



where r^ is the tortoise 

2Mln[r/(2M) - 1]. The quantity V^ OJar/Axial (r) is a po- 
tential that depends on the parity and harmonic number 
I (its precise form can be found elsewhere. To follow our 
notation, see [32j )• 

When the master functions are known, one can use 
them to construct the plus and cross-polarized gravita- 
tional waveforms via 



1 

27 



E 



e>2. 



'(1 + 2)! 
(i-2)\ 



tm 

ZM 



(26) 

where _ 2 Y denotes spherical harmonics of spin weight 
—2 (36[. One can also compute the fluxes of energy and 
angular momentum emitted toward infinity and also into 
the BH horizon in terms of the master functions. These 
fluxes are evaluated using a short-wavelength approxi- 
mation in the radiation zone, where we can introduce 
a well-defined gauge-invariant energy- momentum tensor 
for gravitational radiation [3?], [H, [3!| . However, the ex- 
pression for this effective energy-momentum tensor de- 
pends on the structure of the field equations. In the case 
of CS modified gravity, given that the modification in the 
field equations is additive, the form of these fluxes will 
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be the same as in GR but with extra terms proportional 
to the CS scalar field and its derivatives. That is, 



E, 



1 (1 + 2)! 

64tt ^ (£- 2)1 

O{d6,d 2 0), 



£m 1 2 , |.j-£rn|2 
cpm I 1 r ZM I 



(27) 



L, 



— E 



(1+2)! 

64^ ^ """(£-2)! 

£>2,m v 1 



( \S/ em * 

I CPM 



CPM ^ ZM ZM 



+ 0(d9,d 2 9), (28) 
where the dots here denote time differentiation. 

B. Pontryagin Constraint 

We have seen that diffeomorphism invariance requires 
an extra condition, the Pontryagin constraint [Eqs. ([2])- 
([5])], This condition is automatically satisfied not only 
by the most general spherically symmetric metric (its 
C-tensor vanishes), but also for linear perturbations of 
Minkowski spacetime. Nonetheless, this condition is not 
satisfied for generic perturbations of a Schwarzschild BH. 
In this case, at first-order, the Pontryagin constraint be- 
comes 



*RR 



96M 



•£m 



(h 



r,t 



£(£+ l)Y tm .(29) 



Remarkably, the Pontryagin constraint involves only 
axial modes, which is a consequence of the appear- 
ance of the Levi-Civita tensor (a completely antisym- 
metric tensor) in the modification of the gravitational 
sector of the action. Perhaps even more remarkably, 
this specific combination of axial modes corresponds ex- 
actly to the Cunningham-Price-Moncrief master func- 
tion [331 ] [eq. (|21| ], which appears in the metric wave- 
forms [Eq. (|26p ] and in the fluxes of energy and angular 
momentum [Eq. (|27|) and (|28[)]. We can thus rewrite 
Eq. HUD as 



*RR = - 



2AM (I + 2) 
{£ - 2) 



'_\u£m xr£m 

J ^CPM ■* 



(30) 



Then, the Pontryagin constraint forces the CPM function 
to vanish for all harmonics with i > 2. Such a restriction 
does not imply that all axial perturbations necessarily 
vanish in CS modified gravity, but it does require that 
these modes satisfy the relation 



IL r,t ~ n t,r 



-h e r . 

r 



(31) 



Such a condition seems to reduce the set of possible solu- 
tions of the perturbative vacuum field equations, which 
might lead to an overconstrained system. For the choice 
of 9 in Eq. (p~5|) . we shall see in Sec. IV Al that this is 
indeed the case. 



C. Structure of the Modified Field Equations 

The CS modified dynamics of linear perturbations 
about a Schwarzschild background can be studied by ex- 
panding the metric perturbations into spherical harmon- 
ics, Eqs. (fl"5|) - (fi~7)) . After introducing these expansions 
into the field equations 



C = —C 



(32) 



we can extract individual evolution equations for each 
harmonic. 

The choice of the CS scalar field will determine the 
form of the right-hand side in Eq. (|32[) . Since we are con- 
sidering perturbations about a Schwarzschild spacetime, 
we must choose 9 such that Birkhoff 's theorem holds. In 
Sec. IIII1 we determined that scalar fields of the form of 
Eq. (fT2")l would indeed allow Birkhoff's theorem to per- 
sist in CS modified gravity. We shall here choose 9 as in 
Eq. (TT3")) . such that this field agrees with the symmetries 
of the background, namely 



8(t,r). 



(33) 



where t and r are here the Schwarzschild time and ra- 
dial coordinates. This family of CS coupling functions 
encompasses the canonical choice: 9 = t/fi. 

The structure of the linearized field equations is ana- 
lyzed by first looking at the harmonic decomposition of 
the Einstein tensor and the C-tensor. The structure of 
the former is well-known from the study of perturbations 
of non-rotating BHs in GR, and is given by 



(34) 
(35) 



(36) 

where U^ lar denotes the set of (-£, m)-polar perturbations 

denotes the set of (£, m)-axial perturbations 

(38) 



U AB 


r»lra [rrte 
= VAB l^Polar 


\Y lm , 






^ Aa 


= G'rPZL 




-<wi1i 


gtrn 


u ab 


+ l tm [U l ™ M 


\r£m 
Y ab " 

oirn 
D ab ) 




ryira 
^ab 



and U\ 



l£m 



n__/£rn 

n A , 



Expressions for the coefficients Ga"b> Ga™! G l 
TL lm , and I lm are given in Appendix [A"l Clearly, po- 
lar spherical harmonics have functional coefficients that 
depend only on polar metric perturbations, while axial 
spherical harmonics have functional coefficients that de- 
pend only on axial metric perturbations. On the other 
hand, the harmonic structure of the C-tensor is given by 



t^AB 


p£m \ j jlm 

— ^ABi U Axial 


i Y lm 


film 
^Aa 


nlra \j rim 

— L A [ U Axial 


yim 


r~i£ra 


= C^" l [t/i™al 


y<m 
1 ab 




+ £ em [u e r z ar 


Qirn 
D ab ) 



+ V l A m 



u 



Cm. 



nlm 
D a i 

rytm 
^ab 



(39) 
(40) 

(41) 
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where explicit expressions for the coefficients C A "^, ^a™' 
C em , "D^i T^ lm ') and £ are also given in Appendix [A"l 
In this case, polar spherical harmonics have functional 
coefficients that depend on axial perturbations, while ax- 
ial spherical harmonics have functional coefficients that 
depend on polar perturbations. The main consequence 
of this fact is that in CS modified gravity modes with 
different parity are coupled, and hence, in general they 
cannot be treated separately. 

The linearized field equations, after harmonic decom- 
position, become 



r>tvn 

=J AB = 


nira 
~ U AB J 


Ga = 


nlm 


(42) 


glm _ 


_ c £m ( 




_ v im j 


(43) 


n A = 




j-tm 




(44) 



We can view these equations as the standard Einstein 
equations, linearized about a Schwarzschild background, 
with "source terms" that depend linearly on metric per- 
turbations of opposite parity and their derivatives. Such 
a coupling between different parity modes is analogous to 
what occurs to gravitational- wave perturbations about a 
Minkowski background, where left- and right-polarized 
perturbations mix [J 0, 53] ■ The intrinsic decoupling of 
the modified field equations into polar and axial sectors 
is thus lost. 



and hence the field equations reduce to: 

HT = , l tm = , (48) 

C%=Q, C A m = 0, (49) 

C tm = , V tm = . (50) 

Looking at the expressions of C tm and Cf™ (see Ap- 
pendix [5]) we see that both of them are proportional to 
the metric perturbation h l r m . Therefore, using Eqs. (T4"9"l) - 
(f5TTj) we conclude that 

h e ; n = . (51) 

Similarly, Cj™ is proportional to hl m , and thus, using 
Eq. (T5DD 

h e t m = . (52) 

Since H lm is zero in the Regge- Wheeler gauge, we con- 
clude that all axial perturbations must vanish. In sum- 
mary, the Pontryagin constraint together with the cou- 
pling of opposite parity modes in the modified field equa- 
tions, forbids the existence of purely polar oscillations of 
a Schwarzschild BH in CS modified gravity with a canon- 
ical embedding. 

2. Pure Polar Perturbations 



V. BLACK HOLE PERTURBATION THEORY 
WITH A CANONICAL EMBEDDING 

In order to understand the dynamics of the pertur- 
bations that derives from this theory and, in particular, 
the role and consequences of the Pontryagin constraint, 
we shall concentrate on the special case of a canonical 
CS coupling function, i.e. 6 = t/fi. This canonical cou- 
pling function leads to the canonical timelike embedding 
u„ = [1//^, 0,0,0] and its acceleration ?;„.„ — T c ^ u v a . 



A. One-handed Perturbations 

To begin with, we concentrate on perturbations with a 
single handedness, i.e. purely polar or purely axial per- 
turbations. We shall analyze these cases separately. 



1. Pure Axial Perturbations 

We begin with pure axial metric perturbations, that is 
h% = Pa 1 = K lm = G lm = . (45) 
This conditions imply 



Gab = GT = G £m = H tm = Q, (46) 

V ira _ £ tm = Q ( (47) 



Let us now study pure polar perturbations by setting 
all axial modes to zero: 



h l x i = . 



The immediate consequences are 



Qtm q em Q n 2}^ m Q 



and the field equations become 



nircl nt'm rtlm njim r\ 

y ab = Ga =y = ri =0, 



£ lm = . 

From Ti lm =0we find that (see Appendix [A"|) 



(53) 



(54) 
(55) 



(56) 
(57) 



h l t T = fK? , (58) 

while, from equations Q A m = = Gab we nn d expressions 
for all first and second derivatives of K lm in terms of hf." 1 , 
hf™ and its derivatives. Substituting these expressions 
into £ tm = leads to 



h'™ = , 



which combined with Eq. (|58[) implies 



(59) 



(60) 



There are now only two non-zero metric perturbations: 
K lm and /if™. Inserting the expressions for the deriva- 
tives of K into T> l A l — 0, we can solve the resulting 
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equations for h ™ rr and h t ™ tr . Using all the information 
we have collected so far in Q lm = 0, one can show by di- 
rect evaluation that this equation requires that hr™ t = 0. 
Such a result, in combination with the previously found 
expression for hp? tr , leads to 

K em = 0. (61) 

Returning to the previous expressions for derivatives of 
K lm we obtain 

h*? = 0. (62) 

We have then found that all polar perturbations van- 
ish. We conclude that purely polar oscillations of a 
Schwarzschild BH are not allowed in CS modified gravity 
with a canonical embedding. 

The results obtained for single-parity oscillations are 
quite surprising and raise questions about its robustness. 
In other words, can we expect the same conclusion if 
we repeat the analysis for other CS coupling functions? 
We have repeated this analysis for different CS coupling 
functions within the class 9 = 9(t,r). The algebra in- 
volved is significantly more complicated and it requires 
intensive use of symbolic manipulation software [41] . In 
all studied cases, we have arrived at the same conclusion: 
CS modified gravity does not allow for single-parity BH 
oscillations. 



3. General Perturbations 

General perturbations are significantly more difficult 
to analyze, since we cannot separately study the compo- 
nents of the Einstein and C-tensors, as was the case for 
single parity oscillations, due to the particular structure 
of these tensors [Eqs. (|42]) - (|44l) ]. Nonetheless, we can still 
make some general comments about the consequences of 
the Pontryagin constraint on general oscillations and the 
extent to which these are restricted. 

The Pontryagin constraint unavoidingly adds an extra 
equation to the modified field equations, i.e. Eq. (f3Tj) . 
Since this equation only involves axial perturbations, one 
may a priori think that polar modes are unaffected. The 
modified field equations, however, mix polar and axial 
modes, and thus, polar modes are also affected and re- 
stricted by the Pontryagin constraint. This situation has 
been clearly illustrated by the study of single-parity os- 
cillations in the previous subsection. 

The coupling of polar and axial metric perturbations 
in the modified field equations, together with the extra 
condition imposed by the Pontryagin constraint, lead to 
new conditions on the metric perturbations. In other 
words, these new conditions constitute new equations 
for the metric perturbations that are independent of the 
previous ones. For example, the components of C^S m 
Eq. (|39j) are linear combinations of the axial metric per- 
turbations h g ™ (see Appendix|Aj , which can be combined 
to reconstruct the CPM master function, leading to an 



additional equation for polar modes. Similarly, the com- 
ponents of Tt 1 ™ [Eq. ([35]) ] can also be combined to con- 
struct the CPM master function, and hence find another 
constraint on polar modes. Again, due to the mixing of 
different parity modes, these constraints on polar modes 
can in turn produce new constraints on axial modes. It is 
unclear where this chain of new constraints on the met- 
ric perturbations ends, but it is very likely that they will 
severely restrict the set of allowed BH oscillations. 

One may think that choosing the CS coupling func- 
tion 9 appropriately, namely choosing 9(x ) and Q(x a ) in 
Eq. (|12|) , may end this chain of new constraints through 
cancellations in the generation of new equations. In the 
case 9 = 9(x A ), if such cancellations were to occur for 
a certain harmonic number I, they could not happen for 
other £, as the equations depend on the harmonic num- 
ber. Moreover, such cancellations would have to occur in 
different ways, associated with the different combinations 
of equations that produce new constraints, discussed in 
the previous paragraph. Adding the term <d(x a ) to the 
CS scalar field only complicates the field equations fur- 
ther by inducing a coupling between perturbations with 
different harmonic number. Even in this case, given that 
the functional coefficient of <d(x a ) is fixed (it is just r) 
and ^-independent, it seems unlikely that the coupling of 
harmonic modes will lead to the cancellations necessary 
to avoid new equations for the metric perturbations. 

Our analysis suggests that the present structure of CS 
modified gravity does not seem to allow for generic BH 
oscillations. In particular, we have shown that for certain 
choices of the CS scalar field, single-parity perturbations 
are not allowed. For general perturbations, and with the 
help of computer algebra, we have analyzed the equa- 
tions as discussed above and found that cancellations in 
the generation of new conditions on the metric pertur- 
bations are very unlikely. Therefore, the present set up 
for CS modified gravity does not seem to allow for the 
study of generic oscillations of non-rotating BHs. Such 
a result is reminiscent to that found in GR when addi- 
tional restrictions are imposed on the metric tensor. An 
example of this can be found in the study of relativistic 
cosmological dynamics [H, [43[ 

Although the modified theory seems too restrictive, it 
is possible to develop an extension where the aforemen- 
tioned constraints are avoided, while keeping the main 
characteristics of the modified theory untouched. Two 
distinctive features of the non-extended theory are the 
coupling of different parity modes and the Pontryagin 
constraint. While the former can be preserved by exten- 
sions of the modified theory, the latter must be relaxed 
to prevent the vanishing of the CPM master function. 
Since, by assumption, any extension should introduce 
small modifications to CS gravity, the CPM master func- 
tion would be forced to be small but not quite zero. If 
this is the case, such extensions would have strong im- 
plications in the energy flux and, in particular, in the 
angular momentum flux of gravitational waves [Eqs. (|27|) 
and l|28p]. These arguments motivate the study of BH 
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oscillations in extensions of CS modified gravity, which 
we shall explore in the next section. 



VI. BEYOND THE CANNON 

The strong restrictions obtained on the dynamics of 
BH oscillations thus far can be bypassed provided we 
considered extensions of CS modified gravity beyond the 
canon. One such possibility is to consider a more general 
CS coupling function that, for example, is not spherically 
symmetric. However, as we have argued before, it does 
not seem likely that such a modification would avoid the 
generation of new restrictions on the metric perturba- 
tions. Another more interesting possibility is to extend 
the action of Eq. |T]) to allow for the dynamical evolution 
of the CS scalar field. Such a route is particularly promis- 
ing because it weakens the Pontryagin constraint, as we 
shall see in this section. However, a certain amount of 
arbitrariness is inherent to such a route, encoded in the 
choice of the scalar field action and, in particular, its 
potential. We shall discuss this and other issues in the 
remaining of this section. 



A. Extended CS Modified Gravity 

Until now, we have treated the CS coupling function 
as a non-dynamical quantity. Recently, Smith, et. al. [3| 
added a kinetic and a potential term to the action, which 
they found did not contribute to their weak-field analysis. 
These additional terms are of the form 



1 



V(9) 



(63) 



where V(9) is some potential for the CS scalar field. 

The variation of the extended action with respect to 
the metric and scalar field yields the equations of motion 
of the extended theory 



G>„ + = 87r (T™ at + T, 



[lis ) ' 



no 



dV 
~d9 



R R . 



(64) 
(65) 



where the D'Alambertian of any scalar can be computed 
from 



□0 = 



IsT &A 



(66) 



and where the stress-energy tensor of the scalar field 
is Hi 



rpV 



g^V(9) . (67) 



In this extension of CS modified gravity, the Pontryagin 
constraint is replaced by an equation of motion for the 
CS scalar field, Eq. (|6"5")) . where the quadratic curvature 



scalar *RR plays the role of a driving force. Moreover, 
looking at the field equations for the metric, we realize 
that they are not only sourced by the matter fields but 
also by the CS scalar field through its energy-momentum 
content. In the next subsection we shall study what con- 
sequences this modification imposes on the satisfaction 
of Birkhoff's theorem and the Pontryagin constraint. 



B. Birkhoff's Theorem in the Extended Theory 

The simple derivation of Birkhoff's theorem in Sec. IIIII 
is now modified by the extended action. For a CS scalar 
field of the form of Eq. (| 1 2|) and the line element of 
Eq. (fTQ|) . the C-tensor vanishes, which now leads to 



(68) 



If the metric of Eq. (fT0|) is a solution to Einstein's equa- 
tions, as is the case for the Schwarzschild metric, the 
stress-energy tensor of the scalar field must vanish. For 
the choice of 9 in Eq. (f33|) and the Schwarzschild metric, 
this implies the following conditions: 



= 



° = 2 



n + £V ± fv{9), 



(69) 



where / is given in Eq. (fTTjl . The only solution to this 
system is the trivial one: V = and 9 — const., which 
reduces CS theory to GR. However, if we treat 9 as a 
small quantity, as it is suggested by the different phys- 
ical scenarios that motivate CS modified gravity, then 
Birkhoff's theorem holds to 0(9 2 ) provided V is at least 
of the same order. An example of such a potential is a 
mass term V — m9 2 , typical of scalar interactions. 

The results of this subsection also apply to more 
general CS coupling functions and background metrics. 
In fact, these results hold for any line element that 
represents a general, spherically symmetric spacetime, 
i.e. Eq. (|10p . Moreover, they also hold for the most 
general CS field that leads to a vanishing C-tensor, 

_____ rn\ 

i.e. Eq. (TT2)) , because the kinetic sector of T^J is always 
quadratic in 9. 



C. BH perturbations in the Extended Theory and 
the Pontryagin constraint 

In the extended theory, the Pontryagin constraint is 
replaced by a dynamical equation for the CS scalar field 
with a purely gravitational driving term. Since now there 
is no Pontryagin constraint, there is no a priori reason for 
the equations to disallow general BH oscillations. In fact, 
as we shall see in this subsection, the extended theory 
leads to a system of 11 PDEs for 11 dynamical variables 

But can we treat 9 as a perturbation? As we have just 
seen, Birkhoff's theorem holds only to linear order in 9, 
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for a wide class of potentials V(0). On the one hand, 
in most string-theory scenarios that necessitate the CS 
correction [lO, Q US] , the scalar field 9 is proportional 
to the string scale. In such cases, 9 is much smaller than 
any metric perturbation and BirkhofFs theorem holds. 
Nonetheless, even within such frameworks, the CS correc- 
tion could be enhanced by couplings to non-perturbative, 
string theoretical degrees of freedom (i.e. gravitational 
instantons) (45|. Moreover, there are some theoretical 
frameworks where the strin g coupl ing q s v anishes at late 
times [M El |H, HI M, HHM El Sll , in which case 
a larger coupling is permissible. 

From this discussion, we cannot necessarily assume 
that the magnitude of the scalar field is smaller or of the 
same order as the metric perturbations. In fact, these 
perturbations and the scalar field act on completely in- 
dependent scales. Therefore, all we can assume is that 9 
and | h a b | are both independently smaller than the back- 
ground, which is enough to justify the use of perturbation 
theory. 

We shall thus consider a two-parameter (bivariate) per- 
turbative expansion of CS gravity. One perturbative pa- 
rameter shall be associated with the metric perturba- 
tions, e, and the other with the scalar field, r. The met- 
ric in Eq. (fT4| and the scalar field can then be rewritten 
as 



9 = T(9 + eS9) 

= T9 + re ^ ~9 lm Y lm , 



(70) 
(71) 



e>i,m 



where 9 satisfies BirkhofFs theorem [Eq. (fl"2]) ] and re- 
spects the spherical symmetry of the background. The 
quantities 6 = 9 hn (x A ) are harmonic coefficients of the 
scalar field perturbations, associated with the BH oscil- 
lations. There are no I — modes in the sum of Eq. (|7Tj) 
because they can always be absorbed in the monopolc 
term 9. 

The equations of motion for the metric perturbation 
and the scalar field now become formal bivariate expan- 
sions in e <C 1 as well as t « 1. The modified field 
equations to zeroth order in e are simply equations for 
the background metric, which are automatically satisfied 
to this order by BirkhofFs theorem. The equation of mo- 
tion for the scalar field to the same order becomes 



tU9 = Tg^V M V„fl = . 



(72) 



In Eq. (|72p . and in all remaining equations, the overhead 
bars on any quantity are to remind us that these quanti- 
ties are to be evaluated with respect to the background 
metric g M „. Moreover, in Eq. (|72| and henceforth, we 
shall neglect the contribution of the potential, by assum- 
ing that it is at least of 0(t 2 ). As discussed earlier, the 
potential encodes a certain arbitrariness in the extension 
of the CS modification, which is why we have chosen to 
neglect it. 



The first-order equations govern the dynamics of the 
metric perturbation and the perturbations of the scalar 
field. The equation of motion for the scalar field to 0(e) 
is given by 



e-8(*RR) = e T ^D69- 
- 9„h^„ + 



w + (In y/=g) 



(73) 



In Eq. ([731), 5(*RR) is the functional coefficient of *RR 
to (9(e) [Eq. ([201)] and h = g^V = g^h^. Equa- 
tion ([73]) is harmonically decomposed in the Regge- 
Wheeler gauge in Eq. (|B2|I of Appendix [5] 

Similarly, the equations of motion for the metric per- 
turbation to 0(e) reduce to 



o( 



t 2 ) = eSG 



0,X ap {t N a 5R u 
{V T 9^)8*R J { / V) + {\7 T 59^ 



)0 



(74) 



In Eq. ([74"]) . <5G M „ is the functional coefficient of to 
0(e) (see Appendix |A"]) . and all other 5A a(3 --- x ^,,, stands 
for the coefficient to 0(e) of any tensor A a @--- x <^.... In fact, 
Eq. (|74p is the formal covariant expression of the pertur- 
bation of the Einstein and C-tensors to 0(e), which we 
computed in Appendix [A] for the CS coupling function 
in Eq. (|33[) , with the exception of the term proportional 
to 69. The quadratic terms on the right-hand side of 
Eq. (]74p come from the energy-momentum tensor of the 
scalar field 9 [see Eq. ([M]) ]. When taking the divergence 
of this equation, we recover Eq. ([73]) only when quadratic 
terms in r are taken into account because functional dif- 
ferentiation with respect to 9 reduces the order in r of 
the resulting expression by unity. 

Equations ([73]) and ([74]) govern the dynamics of BH 
oscillations in this extended version of CS modified grav- 
ity, but how do we solve them? Although it would be 
useful to decouple these equations in terms of master 
functions as done in GR, this is not an easy task. In ad- 
dition to the mixing of parity modes and the fact that the 
Pontryagin constraint can no longer be used to simplify 
equations, there is now additional terms in the modified 
field equations due to the perturbations of the scalar field. 
Moreover, these perturbations possess their own dynam- 
ics, and hence, any decoupling should involve the whole 
set of perturbative variables, i.e. (9 im , h e ™). 

Nonetheless, it should be possible to numerically solve 
the perturbative field equations of the extended modified 
theory in an iterative way with r<l. One such possi- 
ble iterative procedure is as follows. To zeroth order in 
r, the equations reduce to those of GR, which can be de- 
coupled and solved numerically with standard methods. 
The numerical result can then be reinserted in the field 
equations to first order in r. These equations can now 
be decoupled in exactly the same way and will contain 
source terms determined by the zeroth order solution. 
This problem will be tackled in a future publication. 
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Before concluding, let us make some general remarks 
about the consequences of extending CS modified grav- 
ity through such kinetic terms. Equation (|73f shows that 
the Pontryagin constraint has become an evolution equa- 
tion for the perturbations of the scalar field, 69. Together 
with Eq. (|74|) . this evolution equation constitutes a sys- 
tem of PDEs for the perturbative variables (8 em ,h e ™) 
that could in principle allow for generic oscillations. How- 
ever, since r -C 1, the magnitude of the CPM master 
function is forced to be small [of 0(t)]. Otherwise, the 
scalar field would lead to an amplification of the CS cor- 
rection to levels forbidden by solar system tests • The 
extended theory thus relaxes the vanishing of the CPM 
function and replaces it by somewhat of a suppression of 
radiative axial modes. 

Such an effect may lead to important observational 
consequences. In particular, the dynamics of astrophysi- 
cal systems, where axial modes contribute significantly to 
the gravitational-wave emission, would be greatly mod- 
ified. In such cases, Eqs. P7|) and (|28p suggest that 
the flux of energy emitted would be dominated by po- 
lar modes and the flux of angular momentum would be 
very small, since it would be linear in r. Such a flux 
suppression is in contrast to predictions of GR, where 
the gravitational wave emission of angular momentum 
is known to be large (approximately 14% of the initial 
ADM angular momentum for quasi-circular BH merg- 
ers [56j]). Consequently, the dynamics of gravitational 
wave sources in the radiation-reaction dominated phase 
should be quite different in CS modified gravity relative 
to GR, thus allowing for gravitational wave tests of the 
extended theory. 



VII. CONCLUSIONS 

We have studied Schwarzschild BH perturbation the- 
ory in CS modified gravity. We began by showing that 
Birkhoff 's theorem (the statement that the Schwarzschild 
solution is the only vacuum, spherically-symmetric so- 
lution of the theory) holds in the modified theory for a 
wide class of CS coupling functions. We then decomposed 
the metric perturbations into tensor spherical harmonics 
and found the linearized modified field equations that 
determine their behavior. The divergence of these equa- 
tions led to the linearized Pontryagin constraint, which 
imposes a restriction on axial metric perturbations -the 
CPM master function has to vanish. 

Once these preliminary issues were studied, we focused 
on the general structure of the metric perturbations in 
CS modified gravity. We found that the modified theory 
adds new terms to the field equations that couple per- 
turbations with polar and axial parity. Moreover, due to 
the restrictions imposed by the Pontryagin constraint, we 
find that in general the system of equations is overdeter- 
mined. Whether the entire set of equations is compatible 
remains unclear, but for a wide class of initial physical 
conditions, we found that linear BH oscillations are not 



allowed in CS modified gravity. Specifically, we showed 
that pure axial or pure polar oscillations are disallowed 
for a wide class of coupling functions. 

Possible extensions of the modified theory that would 
allow for generic BH oscillations were also discussed. In 
particular, we investigated the possibility of providing 
the CS coupling function with dynamics. In other words, 
the inclusion of a kinetic and potential term in the action 
led to the replacement of the Pontryagin constraint with 
an equation of motion for the coupling function. This 
route then lifts the vanishing restriction of the Pontryagin 
constraint and imposes a smallness condition on the CPM 
function. 

The extended CS modified framework thus allows for 
generic BH oscillations but it imposes an important 
smallness restriction on axial perturbation that could 
lead to astrophysical observables. In particular, we saw 
that such a suppression of axial modes would lead to a 
significant decrease in the magnitude of energy, linear and 
angular momentum carried by gravitational radiation rel- 
ative to GR. Such a decrease in gravitational wave inten- 
sity would have important consequences in the dynamics 
of compact object inspirals, specially in the radiation- 
reaction dominated phase. 

Future observations of the ringdown signal in binary 
BH mergers could be used to test and constrain CS mod- 
ified gravity. For such studies, the results found in this 
paper would be critical in order to determine the quasi- 
normal frequency spectrum of perturbations. Future 
work could concentrate on such a spectrum, by numeri- 
cally studying the linearized and harmonically-decoupled 
field equations presented here. Moreover, semi-analytic 
studies might also be possible through the close-limit ap- 
proximation. Only through detailed studies of all aspects 
of the modified theory and its links to experimental ob- 
servations will we be able to determine its viability. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR 
THE PERTURBED C- AND EINSTEIN TENSORS 

Here we provide explicit expressions in Schwarzschild 
coordinates for the coefficients of the expansion in spher- 



ical harmonics of the Einstein and C-tensors. In these 
expressions, the Pontryagin constraint [Eq. 1(9])] has been 
used. For simplicity we omit the superscripts (i,m). 
The components of the harmonically decomposed Ein- 
stein tensor read 



g tt = -fK„ — / ( 7T + — ) K >r 
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while the components of the harmonically decomposed 
C-tensor read 
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APPENDIX B: EQUATIONS OF THE 
EXTENDED CS THEORY IN THE 
REGGE- WHEELER GAUGE 

The equation for the background component of the 
scalar field, Eq. ([72"]). is 

( -figSi + 9 r 2 2 - (r 0) = , (Bl) 
I 



The divergence of the field equations, which in the 
non-extended theory leads to the Pontryagin constraint, 
now leads to the equations of motion for the scalar field. 
To leading order in e, this equation reduces to Eq. (|73| . 
which when harmonically decomposed becomes 
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